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Abstract 

Lecture given at the workshop "Mathematical aspects of theories of 
gravitation" , Stefan Banach International Mathematical Centre, 7 March 
1996. A mini-introduction to critical phenomena in gravitational collapse 
is combined with a more detailed discussion of the regularity of the "crit- 
ical spacetimes" dominating these phenomena. 



1 Critical phenomena in gravitational collapse 

Most of the material here is contained also in the paper M, to which I refer 
the reader for details. Here I'll give detail only in the calculation of the critical 
exponent, and on how gravity regularizes self-similar solutions, where I hope it 
is of more than technical interest. 

Initial data for general relativity (GR) may or may not form a black hole. 
In a hand-waving way one might compare this to a phase transition, and there 
is a "critical surface" in superspace (the phase space of GR) separating the two 
kinds of initial data. Choptuik || has explored this surface in a systematic way. 
For simplicity he took a massless scalar field as his matter model, and allowed 
only spherically symmetric configurations. (The numerical calculations were 
still pioneering work, because they involve length scales spanning many orders 
of magnitude.) To explore the infinite-dimensional phase space, he evolved ini- 
tial data from a number of one-parameter families of data crossing the critical 
surface. Let us call that parameter generically p. By a bisection search, Chop- 
tuik found a critical value p» for each family, such that data with p > form a 
black hole, but not data with p < p*. He then discovered two unforeseen effects: 

Scaling: Near the critical surface, on the black- hole side of it (p > p*), where 
the mass of the black hole final state is small (compared l to, for example, the 
ADM mass), it scales as 

Af ~<7(p-p„) 7 . (1) 
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where the overall factor C depends on the family, but the "critical exponent" 
7 is universal between families. For the scalar field matter, 7 ~ 0.37. I should 
stress that this expression is invariant under redefinitions of p — > p(j>) such that 
p(p) is differentiable, with dp/ dp ^ at p*: To leading order, only C changes 
under such redefinitions. To clarify this, and to avoid speaking in terms of one- 
parameter families altogether, one can formally introduce coordinates (xo,Xj) 
on superspace such that the coordinate surface Xq = is the critical surface. 
Then, for small positive Xo, the black hole mass would be of the form 

M~f(xi)x$, (2) 

and again this form is invariant under coordinate diffeomorphism-s such that 
xo = remains the critical surface. 

Universality: Near the critical surface, on either side of it (p ~ p r ), initial 
data evolve towards an intermediate asymptotic solution (which I'll call Z»), 
which is again universal with respect to initial data. (Clearly this solution 
cannot be a full-blown attractor, because all data must veer off it eventually 
towards forming either a black hole or towards dispersion. But it is an attractor 
of co-dimension one, or intermediate asymptotic.) This solution shows self- 
similarity. In the case of scalar field matter, this self-similarity is discrete, 
showing up as an "echoing" in the logarithm of spacetime scale, with period 
A ~ 3.44. 

The matter model of Choptuik is special in that, even coupled to general 
relativity, it has no intrinsic scale. (This is equivalent to saying that in geometric 
units, c = G = 1, the action has no dimensionful parameters.) Therefore, 
from dimensional analysis alone, there can no be no static, star-like solutions, 
and hence no minimum black hole mass. Generic matter does have both an 
intrinsic scale (or several), and star-like solutions. Astrophysical black holes 
thus probably have a minimum mass given by the Chandrasekhar mass. To 
make infinitesimally small black holes from ordinary matter, one would have to 
use initial conditions (for example rapid implosion) not arising in astrophysics. 
The important point is that for any matter there is at least some region of 
superspace where critical phenomena arise. Nevertheless, the interest of critical 
phenomena does not lie in astrophysics, but in the dynamics of GR (with or 
without matter). Recent work by Choptuik, Bizon and Chmaj ||] (see Piotr 
Bizon's talk, this volume) has clarified the role of matter scales. Investigating 
the spherical collapse of Einstein- Yang-Mills, they find two regimes: one with 
the usual critical phenomena dominated by a self-similar intermediate attractor 
( "second-order phase transition" ) , and one dominated by the Bartnik-McKinnon 
solution, which is a finite mass, static intermediate attractor, and therefore 
having a mass gap ( "first-order phase transition" ) . 

Critical phenomena have also been found in two one-parameter families of 
matter models in spherical symmetry (not to be confused with one-parameter 
families of data!): The first family is that of the perfect fluids with p = kp, 
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k a constant j|, pi k|. The other family is that of the constant-curvature, 
two-dimensional sigma models, characterised by the dimensionless curvature 
parameter re 0. (This family contains the free complex scalar field (§], an 
inflaton-dilaton model |9|, and an axion-dilaton model [[lOj as special cases.) It 
is known for the latter family, and likely for the second, that at some value of 
the parameter re (and perhaps k) the critical solution switches over from con- 
tinuous to discrete self-similarity. Moreover, the critical exponent 7 depends on 
the parameter re or k. 

Historically, the second occurrence of critical phenomena was found in the 
collapse of axisymmctric gravitational waves (time-symmetric initial data, and 
hence with zero angular momentum) fll| . Because of the much greater numeri- 
cal difficulty, there are fewer experimental data, and regrettably, there has been 
no follow-up work so far. Axisymmetric waves are highly interesting in two 
aspects: They go beyond spherical symmetry and they are vacuum data. It 
would also be extremely interesting to see what happens for initial data with 
angular momentum, because black holes can of course have angular momentum, 
but it must be smaller than the mass. 



2 Self-similarity in GR 

The critical solution for each matter model has the property of being self-similar. 
With hindsight one can see that this is the form scale-invariance takes in the 
problem. Self-similarity however takes the novel form of echoing, or discrete 
self-similarity, in some models, notably Choptuik's scalar field and axisymmetric 
gravitational waves. 

The concept of (continuous) self-similarity (CSS) (or homotheticity) has been 
defined in a relativistic context |L2j as the presence of a vector field x sucn that 

C x g ab = 2g ab , (3) 

where C x denotes the Lie derivative. I now introduce the concept of discrete 
self-similarity (DSS). In this symmetry there exist a diffeomorphism <j> and a 
real constant A such that, for any integer n, 

{^) n 9ab = e 2nA g ab , (4) 

where (/>* is the pull-back of <f>. 

To see what DSS looks like in coordinate terms, we introduce coordinates 
(r, x a ), such that if a point p has coordinates (r, x a ), its image 4>(p) has coor- 
dinates (r + A, x a ). One can verify that DSS in these coordinates is equivalent 
to 

g^(T,x a ) = e 2T g^(T, x a ), where <? m „(t, x a ) = g^ v {r + A, x a ) (5) 

In other words, the DSS acts as a discrete isomorphism on the rescaled metric 
g^ v . t is intuitively speaking the logarithm of spacetime scale. 
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In order to clarify the connection between CSS and DSS, one may define a 
vector field \ = d/dr, although there is no unique x associated with a given cf). 
The discrete diffcomorphism cf> is then realized as the Lie dragging along \ by 
a distance A. Clearly, CSS corresponds to DSS for infinitesimally small A, and 
hence for all A, and is in this sense a degenerate case of DSS. In this limit, \ 
becomes unique. 

3 Universality and the critical solution 

The critical solution dominating critical phenomena for a given matter model 
(and perhaps choice of symmetry, such as spherical symmetry) , has two essen- 
tial properties. First, it must be self-similar, either CSS or DSS. Secondly, it 
must have exactly one unstable mode. There are examples of self-similar solu- 
tions with more than one unstable mode, with exactly one (i.e. genuine critical 
solutions), but not with none: the latter would constitute a violation of cosmic 
censorship in the strongest possible sense, that of a naked singularity arising 
from generic initial data. The presence of exactly one unstable mode, on the 
other hand, constitutes a dynamical explanation of the universality of the crit- 
ical exponent. In this explanation, near-critical data p ~ (data near the 
critical surface) are precisely those in which the one unstable mode is initially 
small. Skimming along the critical surface, they are attracted towards the crit- 
ical solution, which is either a fixed point (CSS) or a limit cycle (DSS), until 
the growing mode eventually takes over and ejects the trajectory, either towards 
black hole formation or towards dispersion, in a unique manner. The solution 
has forgotten from which initial data it came, up to one parameter (roughly 
speaking, the number of cchos, or time spent on the intermediate asymptotic) 
which will eventually determine the black hole mass. Figure 1 illustrates this 
behaviour for DSS. 

A given spacetime does not correspond to a unique trajectory in superspace, 
because it can be sliced in different ways. A point in superspace does correspond 
to a unique spacetime, but not to a unique trajectory. In other words, there is no 
preferred time-evolution flow on superspace, and therefore Figure 1 is true only 
for a fixed slicing condition, something extraneous to GR. I am currently trying 
to understand this issue. A possible solution is the fact that the preferred time 
evolution (choice of lapse and shift) should correspond to a "renormalisation 
group flow" , or change of scale. This is the case for example if one takes r 
as defined above as the time variable. The analogy with critical phenomena in 
statistical mechanics would then be deeper: superspace corresponds to the space 
of Hamiltonians, and and the preferred time evolution to the renormalisation 
group flow. 

Constructing the critical solution, one has to proceed in two steps. Impos- 
ing certain regularity conditions and self-similarity, one obtains a non-linear 
eigenvalue problem, which may have a solution. As a second step one has to 
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calculate the spectrum of its linear perturbations and check that there is exactly 
one growing mode. 

In the following I restrict myself to spherical symmetry. The regularity 
conditions just mentioned are imposed in two places. One is the centre of 
spherical symmetry. There one imposes the absence of a conical singularity in 
the 3-metric, and that all fields be either even or odd in r, depending on their 
being vector or scalar under rotations. The other set of regularity conditions 
needs to be discussed in some more detail. 

The general spherically symmetric solution of the wave equation in flat space 

is 

4>(r,t)=r- 1 [f{r-t)+g{r + t)l (6) 

with / and g arbitrary functions. Furthermore, <j) is DSS (in flat space) if 
<p(r, t) = 0(e A r, e A t) for some A. One easily derives that the general DSS 
solution in flat space is of the form 

4>(r,t) = (l + z- 1 ) J F[r + ln(l + z)] + (l-z- 1 )G[r + ln(l-z)], r = ln<, z = r/t, 

(7) 

where F and G are now periodic (with period A), but otherwise arbitrary 
functions. 

With the exception of 4> = 0, a DSS solution can only be regular at either 
r = (F = G), or at r = t (G = 0), or at r = -t (F = 0). All DSS solutions 
(except the zero solution) are singular at the point (r = 0, t = 0). Coupling the 
wave equation to GR in spherical symmetry changes the dynamics, but not the 
degrees of freedom. (There are no spherically symmetric gravitational waves.) 

Surprisingly, the presence of gravity acts as a regulator. In the presence of 
gravity, there is (at least) one DSS solution which is regular at both r = and 
at the past light cone of (r = 0, t = 0) (the generalisation of r = — t to curved 
spacetime) in the sense of being analytic, and regular, in a weaker sense, also at 
the future light cone, (r = 0, t = 0) remains singular, in the sense of a curvature 
singularity. This solution is precisely the intermediate attractor dominating 
near-critical collapse. It is found by numerically solving a non-linear hyperbolic 
eigenvalue problem, with periodicity in some coordinate r, and regularity now 
imposed at both r = and the past light cone. For all details I refer the reader 
to |Tj| |TJ. I'll come back to the regulating effect of gravity in a moment, but 
first I want to conclude my mini-review with an idea of how one calculates the 
critical exponent. 

4 Calculation of the critical exponent 

The calculation of the critical exponent is such a nice piece of dimensional 
analysis [Q, [j| [TJ] that I simply must sketch it here. For simplicity of notation 
I assume CSS of the critical solution, the more generic DSS case is given in the 
paper Q]. Mathematically, the calculation that follows (for CSS) is identical 
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with that of the critical exponent governing the correlation length near the 
critical point in statistical mechanics |D| . 

Let Z stand for the variables of the problem in a first-order formulation, 
in spacetime coordinates and matter variables adapted to the problem jl3| |J. 
Z{r) is an element of the phase space, and Z(r, t) a solution. The self-similar 
solution is of the form Z(r,t) = Z*(r/t). (I won't spell out what Z stands for 
or how r and t are defined.) In the echoing region, where Z* dominates, we 
linearise around it. To linear order, the solution must be of the form 

oo 

Z(r,t) ~ Z* Q + £c < (p)(-t)^ i Z (~) . (8) 

i=l 

Here, the general form of the linear perturbations follows from the form of the 
background solution Z* . Their coefficients d depend in a complicated way on 
the initial data, and hence on p. If Z* is a critical solution, by definition there is 
exactly one Ai with negative real part (in fact it is purely real), say Ai. As t — > 0, 
all other perturbations vanish, and in the following we consider this limit, and 
retain only the perturbation with i = 1. Furthermore, by definition the critical 
solution corresponds to p = p», so we must have C\{p*) = 0. Linearising around 
p*, we obtain 

Jta^O^Z.Q + ^-p.X-t^zg). (9) 

This form holds over a range of t, that is, is an approximate solution. Now 
we extract Cauchy data by picking one particular value of t within that range, 
namely t p defined by 

jr., 

-(p-p*)(-t P ) Xl =e, (10) 



dp 

where e is some constant <C 1, so that at t p the linear approximation is still 
valid. (The suffix p indicates that t p depends on p.) At sufficiently small t, 
the linear perturbation Z\ has grown so that the linear approximation breaks 
down. Later on a black hole forms. The crucial point is that we need not follow 
this evolution in detail. It is sufficient to note that the Cauchy data at t = t p 
depend on r only in the combination r/t p , namely 

Z(r > t p )~Z*(j^)+eS 1 z(j-y (11) 

As furthermore the field equations do not have an intrinsic scale, it follows that 
the solution based on those data must be of the form 

*™ (12) 
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throughout, even when the black hole forms and perturbation theory breaks 
down, and still after it has settled down and the solution no longer depends 
on t. (This solution holds only for t > t p , because in its initial data we have 
neglected the perturbation modes with i > 1, which arc growing, not decaying, 
towards the past.) Because the black hole mass has dimension length, it must 
be proportional to t p , the only length scale in the solution, 



and we have found the critical exponent. 

For completeness I mention that the scaling law is modified when the critical 
solution is DSS. On the straight line relating In M to ln(p — p*) a periodic wiggle 
of small amplitude is superimposed. This wiggle is again universal with respect 
to families of initial data, and there is only one free parameter for each family 
to be adjusted, corresponding to a shift of the wiggly line in the ln(p — p*) 
direction. (No separate shift in the In M direction is required.) 

5 Gravity as a regularizer of self-similar solu- 
tions 

I now come back to the critical solution, which by definition is DSS and regular. 
As I said before, DSS is imposed as periodicity in a coordinate system of the 
form (|5|) adapted to the problem, and boundary conditions at r = arise in 
a straightforward manner, from the necessity of avoiding a conical singularity. 
The generalisation of t = — r to curved spacetime is more interesting. 

From the form (|B|) of the metric it is easy to see that the curvature blows 
up as t -» — oo. Furthermore, this singularity is a "point" in, for example, 
the following sense. Let (n, Ci> #1) ft) an d (t2, C2, #2, V2) be two points. Their 
geodesic distance vanishes as e T as t\ — > T2 — > -co, for any values of (£1, 9±, <p\) 
and (C2,#2,¥>2)- 

The generalisation of r = t to curved spacetime is the past light cone of this 
singularity. We use the remaining freedom in the coordinate system (|^) to label 
this light cone £ = 0. From the analogy with the general solution of the wave 
equation in flat space one would assume that any solution regular at r = is 
singular here, showing an infinite number of oscillations of the form </> ~ i^(ln C) 
(with F periodic) as £ — » 0. 

But this is not so, and I now show why. In flat space, the inward and outward 
travelling modes are cj>' + cf> and </>' — <f>. Let us call their curved-space equivalents 
X + and It is that we expect to be singular at £ = 0. The equation for 
to leading order in £, is 



M oc t p cc {p — p*) A i , 



(13) 



A(t)X_ - X_, T + C(t) 
0B(t) 



(14) 
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where the coefficients A, B and C are constructed from the other fields and are 
therefore periodic in r. 

This approximate equation admits an exact general solution, namely 

x = Jinhom^) + X hom (C, r). (f 5) 

The particular inhomogeneous solution X" lhom is defined as the unique solution 
of 

AX iuhom _ ^inhom + C = Q 

with periodic boundary conditions. This solution exists and is unique, unless 
the average value of A vanishes. The general homogeneous solution X^° m is of 
the form 

x h om = c ^ Ja-^Jb p 

where Aq is the average value of the periodic function A, and J A is its principal 
function after the average value has been subtracted, so that J A is by definition 
also periodic. The periodic function F depends on the initial data for the 
equation (p^). I have not given the expressions for A, B, and C here, but in 
flat space A and B vanish. In the critical solution, which is far from flat, Aq is 
negative and Bo is positive. The exponent of Aq/Bq, is negative for the critical 
or neighbouring solutions. In consequence, X_ m only has two alternatives, it 
either blows up at £ = 0, or it is analytic there (for F = 0). To impose regularity 
at £ = 0, it suffices therefore to impose 

A(t)X_ -AT_, t + C*(t) =0 (18) 

at ( — 0. As it happens, there is locally just one solution which obeys this 
condition, as well as the other boundary conditions. 

As I have said, A vanishes in flat space. But then the equation ( |l8|) has no 
solution X- with periodic boundary conditions, because the average value of C 
does not vanish. But the presence of the term proportional to X^ changes the 
character of the equation quantitatively, and a solution always exists. (In the 
limit as A vanishes, this solution blows up.) 

We have now found the solution in the past light cone of the singularity. The 
data on the past light cone then determine the solution up to the future light 
cone. (We can go no further because the future light cone of the singularity 
is a Cauchy horizon.) Calculating this maximal extension numerically requires 
two more nontrivial changes of coordinate system, giving rise to fresh eigenvalue 
problems. We arrange the final coordinate patch so that once more the future 
light cone is a coordinate line. Now it is X + which is potentially singular. 
Its equation is of the same form (|l4|) , where X + replaces X- , and ( has been 
redefined so that ( — is now the future light cone. Now, however, there is no 
freedom left to adjust any data in order to set F = 0, and in fact F does not 
vanish with the data we have in hand. So the solution cannot be analytic at the 
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Bn 



(17) 
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future light cone. In contrast to the past light cone, however, both Aq and Bq arc 
positive. This means that X^ om , with an infinite number of oscillations as £ = 
is approached, is present, but vanishes at £ = as a (small) positive power of (. 
X + exists at ( — 0, but dX + /d( does not. From the Einstein equations, which 
I have not given here, it follows that the metric and all its first derivatives exist, 
but not some of its second derivatives. Nevertheless these particular second 
derivatives cancel out of all components of the Riemann tensor, so that the 
Riemann tensor exists (but not some of its first derivatives). 

Figure 2 summarises the global situation, with one angular coordinate sup- 
pressed. We have spherically symmetric, discretely self-similar spacetime, with 
a single point-like singularity. The self-similarity corresponds roughly speaking 
to periodicity in the logarithm of the distance from the singularity. The past 
light cone of the singularity is totally regular, in no way distinguished from 
other spacetime points. The future light cone, or Cauchy horizon, is perhaps as 
regular as one can expect, with the scalar matter field C°, the metric C 1 and 
the Riemann tensor C . In particular it carries well-defined null data (which 
are of course self-similar), and there exists a regular, self-similar (and of course, 
non-unique) extension of the spacetime inside the Cauchy horizon, with only 
the horizon itself of limited differentiability. Surprisingly, the null data on the 
horizon are very small, so that the extension can be made almost flat and empty. 
The situation is similar in the other two cases where the critical solution has 
been calculated up to the Cauchy horizon [|[ [H^. The spacetime there is CSS, 
but also almost flat at the horizon. A limiting case arises in the closed-form 
solution of Roberts fHI , where the null data on the horizon vanish exactly, 
so that a possible extension inside the light cone is flat empty space. (In the 
Roberts solution, the past light cone also carries zero data, in contrast to the 
other examples, where the past light cone is very far from flat.) Finally I should 
add that at large spacelike distances from the singularity, spacetime becomes 
asymptotically conical, with a constant defect solid angle. 

Of course, all these global considerations are not relevant to critical collapse, 
where the spacetime asymptotes the critical spacetime in some bounded region 
inside the past light cone. Therefore, neither the singularity nor the asymptot- 
ically conical region, nor the Cauchy horizon appear. 

6 Conclusions 

Looking back on what has become a small industry in the past two years, I 
think that the dynamical mechanism of universality, scaling and echoing is now 
understood in an intuitive way. Furthermore we can calculate the echoing period 
A and critical exponent 7 as nonlinear eigenvalues, in a manner distinct from 
that of fine-tuning data in numerical experiment. 

Critical phenomena have also given new material to the study of cosmic 
censorship. There are self-similar spacetimes (spherically symmetric, except for 
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one axisymmetric example) with a naked singularity that have an infinity of de- 
caying linear perturbation modes opposed to only one increasing perturbation 
mode. This still implies cosmic censorship in the sense that a generic perturba- 
tion, which will contain some small fraction of that unstable mode, destroys the 
naked singularity (either by forming a horizon or by avoiding a singularity alto- 
gether, depending on the sign of the mode amplitude). But we are very close to 
a violation, in that we only have to set one mode out of an infinity of modes to 
zero to get the naked singularity. A single generic parameter in the initial data 
provides a sufficient handle to do this. (In the notation above, we only have to 
set the one amplitude C\ equal to zero, and if C\ depends in some way on p, we 
can do this by adjusting p.) Therefore by arbitrary fine-tuning of one generic 
parameter in the initial data one can obtain asymptotically flat spacetimes in 
which a region of arbitrarily large curvature is visible to an observer at infinity. 

There are a number of important open questions. Are there additional un- 
stable perturbations of the critical solution among the nonspherical modes? Are 
there critical solutions in axisymmetry, or even lower symmetry? Is DSS more 
generic than CSS? How do the charge and angular momentum of the black hole 
scale as one fine-tunes initial data with charge and angular momentum, with the 
aim of making a black hole of infinitesimal mass (and not caring about charge 
and angular momentum)? Is there a preferred flow on superspace which would 
complete the analogy with the renormalisation group flow in statistical mechan- 
ics? Is there even statistical physics hidden somewhere? (Very, very unlikely, 
but who would have thought it of black holes pre-1974?) Can someone give an 
existence proof for the critical solutions, which so far have only been constructed 
numerically? Can one prove that they must have the limited differentiability 
(metric C 1 etc.) at the Cauchy horizon which I have described, not more or 
less? 
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Figure 1: Figure 1: The phase space picture for discrete self-similarity. The 
plane represents the critical surface. (In reality this is a hypersurface of co- 
dimension one in an infinite-dimensional space.) The circle (fat unbroken line) 
is the limit cycle representing the critical solution. The thin unbroken curves 
are spacetimes attracted to it. The dashed curves are spacetimes repelled from 
it. There are two families of such curves, labeled by one periodic parameter, 
one forming a black hole, the other dispersing to infinity. Only one member of 
each family is shown. 
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